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Processing Seismic Data 

The present invention relates to a method of processing multi-component seismic data. 
It particularly relates to a method of processing seismic data to determine a calibration 
filter that calibrates one component of the seismic data relative to another component of 
the seismic data. The invration further relates to an apparatus for processing seismic 
data. 

Figure 1 is a schematic view of a seismic surveying arrangement. In this figure the 
surveying arrangement is a marine surveying arrangement in which seisnaic energy is 
emitted by a seismic source 1 that is suspended within a water column 2 from a towing 
vessel 3. In this exanople the water column is the sea, but the methods described 
hereinbelow can be applied to data acquired in seawater or in freshwater. When the 
seismic source 1 is actuated seismic energy is emitted downwards and is detected by an 
array of seismic receivers 4 disposed on the seafioor 5. (As used herein the term 
"seabed" denotes the earth's interior, and the term "seafloof denotes the surface of the 
seabed.) 

Many seismic surveys now use multi-component receivers that record two or more 
components of the seismic energy incident on the receiver. For example a 3-component 
(3-C) seismic receiver contains three orthogonal geophones and so can record the x-, y- 
and z-components of the particle motion at the receiver (the particle motion may be the 
particle displacement, particle velocity or particle acceleration or even, in principle, a 
higher derivative of the particle displacement). In a marine seismic survey a 4- 
component (4-C) seismic receiver can alternatively be used. A 4-C receiver contains a 
(dynanodc) pressme sensor such as a hydrophone in addition to three orthogonal 
geophones and so can record pressure fluctuations as acoustic waves propagate in the 
water colunm (a scalar quantity) in addition to tiie x-, y- and z-components of the 
particle motion of the seabed. 
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Many different paths exist by which seismic energy may travel from the source 1 to a 
receiver 4 in the seismic surveying arrangement of Figure 1. A nunober of paths are 
indicated schematically in Figure 1. 

The path 6 shown in Figure 1 is known as the "direct path". Seismic energy that travels 
along the direct path 6 travels from the source 1 to a receiver 4 essentially in a straight 
line without undergoing reflection at any interface. 

Path 7 in Figure 1 is an example of a ^'water layer multiple path". Seismic energy that 
follows a water layer multiple path propagates wholly within the water column 2, but 
undergoes one or more reflections at the surface of the water column and/or the seafloor 
5 so that the seismic energy passes tJirough the water column more than once. The 
water layer multiple path 7 shown in Figure 1 involves one reflection at the seafloor 5 
and one reflection at the surface of the water column, but many other water layer 
multiple paths exist. 

The path 8 in Figure 1 is an example of a "critical refraction path". Seismic energy that 
follows the path 8 propagates downwards to the seafloor 5, and penetrates into the 
earth's interior 10 (ie into the seabed). The seismic energy continues propagating 
downwardly, until it reaches a boundary 11 between two layers of the earth that have 
different acoustic impedance. The seismic energy undergoes critical refraction, 
propagates along the boimdary 11, before eventually being refracted upwards towards 
the receiver 4. Critical refraction may also occur at the water-seabed interface, and 
downwardly propagating seismic energy that is refracted in this way wiU propagate 
along the water-seabed interface and will then propagate upwardly into the water 
colmnn. 

The path 9 shown in Figure 1 is known as a "primary reflection path". Seismic energy 
that follows the primary reflection path 9 propagates downwards througih the water 
column, is refracted at the seafloor 5, and propagates downwardly through the earth's 
interior. The seismic energy is refracted at the boundary 11, but is not critically 
refracted and so continues to propagate downwardly into the earth. It eventually 
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undergoes reflection at a geological structure 12 that acts as a partial reflector of seismic 
energy, and the reflected seismic energy is, after further refraction as it passes upwardly 
through the boundary 11, incident on the receiver 4. The general intent of a seismic 
survey is to make use of the seismic energy that follows the primary reflection path in 
order to obtain information about the interior structure of the earth. 

Seisnoic energy acquired at a receiver may contain upwardly and/or downwardly 
propagating seismic energy depending on the location of the receiver and on the event. 
For example seismic energy that travels along the critical refraction path 8 shown in 
Hgure 1 win, when it is incident (travelling upwardly) on the water-seabed interface, be 
partly transmitted into the water colunm and partially reflected back into the seabed 10. 
Thus, a critical refraction event will consist purely of upwardly propagating seisnoic 
energy above the seafloor 5, but will contain both upwardly and downwardly 
propagating seisnaic energy below the seafloor 5, As another example, seismic energy 
that travels along the direct path 6 shown in Figure 1 will, when incident on the water- 
seabed interface 5, be partially transmitted into the seabed and partially reflected back 
into the water column. Hence, the direct event will contain both upwardly and 
downwardly propagating seismic energy above the seafloor, but will contain only 
downwardly propagating seismic energy below the seafloor. It is therefore often of 
interest to decompose the seismic data acquired at the receiver 4 into an up-going 
constituent and a down-going constituent, above or below the seafloor 5. For example, 
in a 4-C seismic survey it may be of interest to decompose the pressure and the vertical 
particle velocity recorded at the receiver into their up-going and down-going 
• constituents above the seafloor. 

Various filters that enable decomposition of seismic data into up-going and down-going 
constituents have been proposed. For example, K,M. Schalkwijk et al have suggested, 
in "Application of Two-Step Decomposition to Multi-Component Ocean-Bottom Data: 
Theory and Case Stud/', J, Seism. Expl. Vol. 8 pp261-278 (1999), that the down-going 
and up-going constituents of the pressure jxist above the seafloor may be expressed as: 
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(1) 



where P is the pressure acquired at the receiver, P' is the up-going constituent of the 
pressure above the seafloor, is the down-going constituent of the pressure above the 
seafloor, / is the frequency, k is the horizontal wavenumber, is the vertical particle 
velocity component acquired at the receiver, p is the density of the water, and q is the 
vertical slowness in the water layer. 

As can be seen, the expressions in equation (1) require two of the conaponents of 
seismic data recorded at the receiver to be combined. These jOlters are an example 
wheie it is necessary to combine two components of the acquired seisnoic data. It may 
also be necessary to combine two or more components of the acquired seismic data in 
order to decompose the acquired seismic data into p-wave and s-wave (pressure-wave 
and shear-wave) components, or to remove water level multiple events from the seisnoic 
data. 

One problem in combining different components of the seismic data acquired at a 
receiver is that the different components of the seismic data may not be correctly 
calibrated against one another. IMs is particularly the case where the two components 
that are being combined are, as in equation (1), the pressure and the vertical particle 
velocity. There are usually differences in coupling or impulse response between the 
hydrophone used to acquire the pressure and the geophone used to acquire the vertical 
particle velocity. It is necessary to calibrate the data for these differences before the 
pressure and vertical particle velocity can be combined, and this process is known as 
'T/vz calibration". This calibration process involves developing a calibration jBlter that 
compensates for the differences in coupling and impulse response between the 
hydrophone and the vertical geophone and then applying the filter to one data set to 
compensate for the differences in coupling. 
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Schalkwijk et al, and others, have suggested that the P/vz calibration problem can be 
addressed by assuming that one component of the seisnaic data has been correctly 
recorded, and calibrating the other component of the seismic data against the component 
that is assumed to be correctly recorded In general, it is assxamed that the hydrophone 
is well coupled to the seismic wavefield, so that the pressure recording is taken to be 
correct. The vertical component of the particle velocity is then calibrated against the 
pressure to compensate for coupling and impulse response differences bet\veen the 
hydrophone and the vertical geophone. Schalkwijk et al therefore proposed that 
equation (1) above should be modified by applying a calibration filter to the vertical 
particle velocity. They proposed that the equation given above for the down-going 
constituent of the pressure above the seafloor should be modified to read as follows: 

i'*(/.«4p(/.*)+<.(/)^^v.(/.*). (2) 

In equation (2) a(f) represents a frequency-dependent calibration filter. The remaining 
terms in equation (2) have the same meaning as in equation (1). 

The method proposed by Schalkwijk et al. for determining the calibration filter is to 
find the calibration filter that minimises the energy of the down-going pressure 
constituent above the seafloor for a pordon of the seismic data that contains only 
primary reflections. Seismic energy travelling along a primary reflection path is 
propagating upwardly just above the seafloor at the receiver position, so that the down- 
going constituent of the pressure just above the seafloor should be zero for a primary 
reflections. Schalkwijk et al. used a least squares method to find the calibration filter 
that minimises the energy of the down-going constituent of the pressure in a window 
containing only primary reflection events. Once the calibration filter a(f) has been 
determined in tibis way, it is applied to the entire data set to calibrate the vertical particle 
velocity. 

The existence of various paths of seisnaic energy from the source to the receiver means 
that the data acquired at the receiver in a real seismic survey will contain events 
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coiresponding to more than one possible path. These events will occur at different 
times after the actuation of the seismic source 1, as different paths of seismic energy 
have different associated travel times. In order to apply the method of Schalkwijk et al. 
to determine the calibration jBlter, data in a time window that contains only a primary 
reflection event must be selected, 

The method proposed by Schallcwijk et al. has the disadvantage that the time window 
containing only primary reflection events has to be picked manually. The primary 
reflection events are not the £b:st events acquired at the receiver following actuation of 
the source, and so cannot be picked automatically. A further disadvantage is that in 
some cases, for example if the seismic source has a long signature, it may be hard to 
distinguish between the direct arrival and the primary reflection events, so that it may be 
difficult to isolate the correct events. Moreover, in shallow water the water layer 
multiple events may arrive shortly after the direct wave. In this case, the derivation of 
a(f) is based on a very limited amount of data, reducing the accuracy of the results. The 
direct event and water multiple events contain downwardly propagating seismic energy 
so that use of a time window that inadvertently included the direct event or water 
multiple events would not give correct results for the calibration filter. 

Co-pending UK Patent application No. 0200560.1 and PCT application PCT/GB 
03/00052 propose applying the approach of Schalkwijk et al to a time window 
containing critically refracted waves, which also contain only up-going energy above 
the seafioor. This technique is more suitable for automation using first-break pickers, 
and allows the method to be applied to data acquired in shallow waters. 

BaU, V. L. and Corrigan, D, suggest, in *T)ual sensor summation of noisy ocean-bottom 
data'*, 66th Ann. Intemat. Mtg: Soc. of Expl. Geophys., 28-31 (1996), caUbrating the 
vertical geophone data against pressure by applying a hydrophone ghost operator to Vz 
and a geophone ^ost operator to flie pressure. It can be shown that this approach is 
equivalent to predicting the down-going pressure reflected from the sea surface from the 
computed up-going pressure. Miniznising the difference between the predicted and the 
computed down-going pressure reflected from the sea surface at times larger than 
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source duration plus the one-way propagation time through the water layer then allows 
the desired calibration filter to be determined. A significant drawback of this techmque 
is the need for accurate information about the depth of the water colunm in the survey 
area. 

Another approach to calibrating data for differences before the pressure and vertical 
particle velocity is to require seismic energy to be preserved during propagation through 
the water layer. In this case, frequency- and wavenumber-dependent calibration 
operators are designed by means of spectral balancing of the up- and down-going wave 
constituents just above the seafloor. 

The present invention provides a method of processing multi-component seisnoic data 
acquired at a receiver station from seismic signals propagating in a medium, the method 
comprising tiie steps of: selecting a first portion of the seismic data; and determining a 
calibration filter firom the first portion of the seismic data, the calibration filter being to 
cahbrate a first component of the seismic data relative to a second component of the 
seisnoic data; wherein the step of determining the calibration filter comprises processing 
the data in the conmion shot domain. 

Hie term "receiver station" may denote, for example, a single receiver (for example in 
singile sensor seismic acquisition), a hard-wired group of receivers, etc. 

Coupling and instrument response variations for a particular seismic receiver or receiver 
station are receivCT-consistrat effects - tiiat is, once a multi-component receiver has 
been deployed the required calibration operator for the receiver is constant and does not 
depend on the source position. Receiver calibration has therefore hitherto been 
performed in the common-receiver domain (CRD). When seismic data is processed ia 
the common-receiver domain, all data acquired at one receiver is sorted into a "gather" 
for that receiver - and since all data in the receiver gather was acquired at a single 
receiver, a single calibration operator can be used to correct all the data in the gather for 
coupling and instrument response variations. 
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The inventors have realised, however, that optimising the vector fidelity of multi- 
component receivers in the common shot domain (CSD) has considerable advantages, 
including: 

unlike CRD processing, processing in the CSD does not implicitly assume the 
sub-surface to be laterally invariant; 

rough sea perturbations can properly be accounted for; 

it is possible to benefit from the dense receiver spacing if the source-side 
spacing is coarser; and 

waterbome noise that is not shot-generated can be removed more efficiently. 

It is conventional to perform tiie calibration technique in the coromon-receiver domain 
since this allows the calibration filter for each receiver to be determined separately, but 
this has the disadvantage that rough sea perturbations cannot be properly accounted for, 
and that background noise cannot be removed efficientiy. Furthermore, applying the 
processing in the common receiver domain implicitiy assumes that the earth's sub- 
surface is laterally invariant 

The calibration filter may then applied to seismic data recorded at the specific receiver 
or receiver station which acquired the seismic data used to determine the calibration 
filter. Applying the calibration filter compensates for the effects of the different 
coupling for the first and second components of the seisnoic data. 

The calibration filter is receiver-specific, so a separate calibration filter is preferably 
detemained for each receiver (or receiver station) in a seismic survey. 

A second aspect of the invention provides a method of seismic surveying comprisiag 
the steps of: actuating a source of seismic energy; acquiring seismic data at a receiver 
station spatially separated from the source; and processing the data by a method as 
defined in the first aspect of the invention. 
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A third aspect of the invention provides an apparatus for processing multi-component 
seisioic data acquired at a receiver station from seismic signals propagating in a 
medium, the apparatus comprising: means for determining a caKbration filter from a 
first portion of the seismic data, the calibration filter being to calibrate a first component 
of the seismic data relative to a second component of the seismic data; wherein the 
apparatus is adapted to determine the calibration filter by processing the seismic data in 
the common shot domain. 

Hie apparatus may comprise a programmable data processor. 

A fourth aspect of the invention provides a storage mediimi containing a program for 
controlling a progranunable data processor to carry out a method as defiined in the first 
or second aspect. 

A fifth aspect of the invention provides a program for controlling a computer to cany 
out a method as defined in the first or second aspect. 

Preferred features of the invention are set out in the dependent claims. 

Preferred embodiments of the present invention will now be described by way of 
illustrative example with reference to the accompanying figures in which: 

Figure 1 is a schematic illustration of a seismic survey; 

Figure 2 is a block flow diagram of a method according to one embodiment of the 
present invention; 

Figure 3 is a block flow diagram of a method according to a second embodiment of the 
present invention; and 

Figure 4 is a schematic block diagram of an apparatus according to the present 
invention. 
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As outliaed above, existing P/Vz calibration techniques are based on finding calibration 
filters that niimmise the up-going or down-going constituent of pressure in a water 
colimm (as defined by equation (2)). The term l/q(f,k) in equation (2) may be 
considered as a spatial filter, and the prior methods have the disadvantage that they 
require computing the term l/q(f,k) [a(f) y^,k)'\. As the calibration filters are a function 
of position, spatial filter operations on the term [a(f) V2(f JO] complicate the 
computation of the calibration filters a(f) for each individual receiver position. 

The present invention overcomes these disadvantages by perfomiing the spatial filter 
operation in the common shot domain. This allows the optimisation criteria to be 
applied to the up- and down-going components of particle velocity, thus avoiding 
spatial filter operations on the term [a(f) vz(f,k)]. 

The invention makes use of the fact that the up- and down-going constituents of the 
vertical component of particle velocity in the water layer can be expressed as 



In equation (3), aJ(f} again denotes the frequency-dependent calibration filter that 
corrects for imperfections in the recording of v^, Vz" denotes the up-going constituent of 
the vertical component of particle velocity, and Vx*" denotes the down-going constituent 
of the vertical component of particle velocity just above the seafloor. The remaining 
quantities in equation (3) have the same meaning as in equation (2). 

The determination of the up- and down-going constituents of the vertical component of 
the particle velocity involves applying the spatial filter q(f,k) to the pressure, rather than 
to the calibrated vertical particle velocity. As a result, the filter operation does not 
interact with the unknown calibration filter aj(f) when applied in equation (3) in the 
common shot domain. Determining the calibration filter in the common shot domain — 
or performing at least some of the steps of the determination of the calibration filter in 




2p 



P(f.k) 



(3) 
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the coiomon shot domain - allows the calibration filters to be determined from the up- 
and down-going components of the particle velocity according to equation (3), therefore 
greatly simplifying the processing. 

A further advantage of using equation (3) is that the filter operator l/q(f,k) used in 
equation (2) contains a pole. In contrast, the filter operator q(fik) in equation (3) does 
not contain a pole, and simply contains a zero. As a result, techniques based on 
equation (3) are numerically more stable, and facilitate the use of spatially compact 
filter approximations. The use of spatially compact filter approximations is 
advantageous when processing data acquired at a survey location where there are 
significant lateral variations in the earth's surface parameters. 

Figure 2 illustrates one method of detemiining a calibration filter in the coiomon shot 
domain. This example is described with reference to P/Vz calibration, but a similar 
procedure may also be carried out for P/Vx calibration as suggested by Schalkwijk et al, 
(above), provided that the elastic properties of the seabed are known. 

The method shown in figure 2 is based on the P/Vz calibration techniques proposed by 
Schalkwijk et al (above) and in UK patent application No. 0200560.1, In the method of 
figure 2, the firequency-dependrat calibration filter is determined by miniroising the 
energy of the down-going constituent of the vertical component of the particle velocity 
Vz" over a data window that contains only up-going events. For example, the data 
window may contain only the primary reflection event (as proposed by Schalkwijk et 
al.) or it may contain only critically refiracted events (as proposed in UK patent 
application No. 0200560.1). 

Initially, multi-component seismic data are acquired at step 20. In this embodiment the 
multi-component seismic data contain at least pressure recordings and vertical particle 
velocity recordings. 

The invention may alternatively be applied to pre-existing seismic data. In this case, 
step 20 is replaced by the step, step 21, of retrieving multi-component seismic data 
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contaiiung at least pressure recordings and vertical particle velocity recordings firom 
storage. 

Next, at step 22, the pressure data is sorted into common shot gathers. Each gather will 
contain all pressure data recorded for a particular shot - so where data is acquired using 
the seismic surveying arrangement shown in figure 1, for example, one gather will 
contain pressure data recorded when the source 1 is actuated at one location indicated 
schematically by 13 in Figure 1, another gather will contain all pressure data recorded 
when the source is at another location indicated schematically by 13 a and so on. Each 
gather will contain pressiire data acquired at all receivers 4 in the seismic surveying 
arrangement. The pressure data acquired by a receiver will typically be in the form of a 
record of the pressure at a receiver measured as a function of time since the actuation of 
the seismic source 1. That is, each receiver records P(receiver position, t) which, if the 
sources and receivers are arranged along a straight line, may be simplified to P(x,t). 

Next, at step 23, the pressure data are transformed, in this embodiment into the 
firequency-wave number (f-k) domain to give P(f,k). The transform to the frequency- 
wave number domain may be carried out using any suitable technique. The invention is 
not, however limited to the f-k domain and the method of the invention may be used in 
other domains such as, for example, the x-t domain, the x-/ domain, the k-t domain, the 
A:-/domain, the tau-p domain, etc. 

Next, at step 24, the pressure data are spatially filtered. That is, the filter q(fik) is 
applied to the pressure data to yield the filtered pressure data q(f,k)P(f,k). The filter 
q(f,k) is obtained firom the seismic data acquired at step 20 or retrieved at step 21, and 
this can be done using any known technique. The determination of the filtered pressure 
data may also be calculated using any known technique. This spatial filtering step is 
carried out on each pressure recording in a shot gather, and this process is then repeated 
for all shot gathers. 

At step 25, the filtered pressure recordings are transformed back to the x-domain, for 
example the f-x, t-x or t-x domain. The filtered pressure recordings are then sorted into 
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receiver gathers. Each receiver gather contains (filtered) pressured recordings acquired 
at only one receiver or receiver station. 

At step 26, a frequency dependent calibration filter a(f) is determined for one receiver 
gather from a portion of the seismic data in that gather. In this embodiment the 
caUbration fUter is determined using equation (3), by finding the calibration filter that 
minimises the up-going or down-going constituent of the vertical particle velocity for a 
time window in which that constituent is expected to be zraro. (Although equation (3) is 
written in the f-k domain, a corresponding equation exists in the x-domain.) For 
example, the time window selected may be a time window that is expected to contain 
only up-going seismic energy just above the sea floor, such as a time window that 
contains only the primary reflection events or a time window that contains only 
critically refracted events. Where a time window that is expected to contain only up- 
going seismic energy is selected, the calibration filter is determined by finding tiie filter 
that minimises the down-going constituent of the vertical particle velodty in tiie 
selected time window. 

The result of step 26 is a calibration filter for tiie particular receiver (or receiver station) 
used to acquire the pressure data in tiie receiver gatiier on which step 26 was performed. 
This calibration filter may tiien be used to calibrate vertical velocity component data 
acquired at that receiver (or receiver station), at step 27. The calibration filter 
determined for a receiver (or receiver station) may be applied to vertical velocity 
component data outside the time window used to determine the calibration filter, and it 
may be used to calibrate vertical velocity component data which was acquired at that 
receiver (or receiver station) but which was not used in the dettamination of the 
calibration filta:. 

Steps 26 and 27 may then be repeated for other receiver gatiiers, to determine 
calibration filters for each receiver (or receiver station) in tiie seismic surveying 
arrangement. 
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The caKbrated vertical velocity data may be subjected to further processing steps (not 
shown). 

Ici the embodiment of figure 2, step 22 comprises transforming the pressure data to the 
f-k domain. It would alternatively be possible for step 22 to comprise transforming the 
pressure data to the x-p domain, in which case step 23 would comprise filtering the 
pressure data in the T-p domain. 

A second embodiment of the present invention is illustrated in figure 3 of the present 
application. This embodiment is based on the principle that the water column is non- 
attenuative for seismic waves, and that the free surface of the water column has a 
known, constant reflection coefficient ro. As a result, all up-going energy in the water 
column should, at some later time, be recorded as down-going energy but with an 
amplitude reduced by a factor ro. The only exception to this principle is the direct wave 
6 - the direct wave 6 produces down-going energy in the water colurmi that does not 
have any corresponding up-going energy. A receiver-consistent, frequency-dependent 
calibration filter can be determined by finding the calibration filter that provides spectral 
balancing between the up-going energy and the down-going energy (after removing, or 
"muting", the direct wave from the down-going energy). Thus, the calibration filter is 
determined as the calibration filter that noinimises the following objective function: 



]y;{f^k))\- \vr''{f^k))\ )\ (4) 

In equation (4), W(f,k) is a weighting function that specifies the frequency-wave number 
window over which the miimnisation process is performed, Vz is the up-going 
constituent of the vertical particle velocity, and y/'"^'*' is tiie down-going constituent of 
the vertical particle velocity after muting of the direct wave. 

The terms Vz(f,k) and V^^^^ (f,k) in equation (4) may be expanded using equation (3) to 
give: 
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2p 2 2p z J 



The superscript (r) again indicates that the direct wave has been removed. 

Muting the direct wave is a known, straightforward process. In essence, a data mask is 
constructed that has a value of zero for the approximate arrival times of the direct wave 
and that is equal to one at all other times. The direct wave is then muted by 
multiplication of this mask with the data. The arrival time of the direct wave can be 
estimated ftom a rough estimate of the survey geometry, water depth and the velocity of 
seismic energy. The transition between zero and one in the data mask is preferably 
smooth rather than abmpt, and can be anywhere between the arrival time of the direct 
wave itself and the arrival time of the first water layer multiple since any primary 
reflections that may occur between these times are purely up-going waves. The direct 
wave is most easUy muted in tiie (x-t) or (tau-p)-domains - in which case the pressure 
data and vertical velocity data are muted and flien transformed to ih&f-k domain. 

Figure 3 is a flow chart showing one example of this method. Initially, multi- 
component seismic data are acquired at step 30 or are retiieved from storage at step 31. 
The multi-component seismic data contain at least pressure recordings and vertical 
particle velocity recordings. 

At step 32, the pressure recordings are sorted into conmion shot gathers. At step 33 the 
pressure recordings are tiansformed to tiie f-k domain, and at step 34 the filter q(f.k) is 
obtained and the pressure data ate filtered to detennine the quantity q(f.k)P(f,k). Steps 
30-34 correspond generally to steps 20-24 of the method of figure 2. 



At step 35, the filtered pressure data are sorted into common receiver gathers. 
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One particular receiver gather is then selected at step 36. The pressure data P(x,t) and 
the vertical velocity data v^x,t) for that gather are then muted to remove the direct 
wave. The muted pressure data and vertical velocity data are then transformed to the f-k 
domain, at step 37, to give P^''^ (f,k) and y/^^ (f,k). The un-muted vertical velocity data 
are also transformed to the f-k domain in this step, to give Vz (f,k). 

At step 38, a frequency-dependent calibration filter is determined for the selected 
receiver gather. This calibration filter is determined as the calibration filter that 
mininoises the objective function E of equation (4) or equation (S). 

The calibration filter obtained in step 38 may tiien be used to calibrate vertical velocity 
data acquired by the receiver (or receiver station) corresponding to the selected gather. 
This is indicated as step 39 in figure 3. 

Steps 36-39 are then repeated for other receiver gathers, to determine calibration filters 
for each receiver. 

The calibrated vertical velocity data may be subjected to further processing steps (not 
shown). 

Figures 2 and 3 illustrate a calibration filter being obtained for one receiver gather (step 
26, 38), and this calibration filter is applied to data for that receiver (step 27, 39) before 
calibration filters are determined for other receiver gathers. The invention need not be 
performed in this way, and it would be possible to determine calibration filters for each 
receiver gather before applying any of the calibration filters to data acquired by the 
respective receiver. 

The invention has been described above with reference to calibrating the vertical 
velocity component Vz against pressure, but the invention is not limited to this. It may 
alternatively be used to calibrate a horizontal component of the velocity against the 
pressure. 
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Hgute 4 is a schematic block diagram of a progjcammable apparatus 14 according to the 
present invention. Hie apparatus comprises a programmable data processor 15 with a 
programme memory 16, for instance in the form of a read-only memory (ROM), storing 
a programme for controlling the data processor 15 to perform any of the processing 
methods described above. The apparatus further comprises non-volatile read/write 
memory 17 for storing, for example, any data which must be retained in the absence of 
power supply. A "working" or scratch pad memory for the data processor is provided 
by a random access memory (RAM) 18. An input interface 19 is provided, for instance 
for receiving conmiands and data. An output interface 20 is provided, for instance for 
displaying information relating to the progress and result of the method. Seismic data 
for processing may be supplied via the input mterface 19, or may alternatively be 
retrieved £rom a machine-readable data store 21. 

The programme for operating the system and for performing tiie method described 
hereinbefore is stored in the programme memory 16, which may be embodied as a semi- 
conductor memory, for instance of the well-known ROM type. However, the 
programme may be stored in any other suitable storage medium, such as magnetic data 
carrier 16a, such as a "floppy disk" or CD-ROM 16b. 



